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Abstract 

We investigate two kind of polarization of localized optical waves in 
nonlinear Kerr type media, linear and combination of linear and circu- 
lar. In the first case of linear polarized components we have obtained 
the vector version of 3D+1 Nonlinear Schrodinger Equation (VNSE). 
We show that these equations admit exact vortex solutions with spin 
£ = 1. We have determined the dispersion region and medium param- 
eters necessary for experimental observation of these vortices in the 
conlusion. In the second case we represent the electric and magnetic 
fields as a sum of circular and linear components. We supouse also 
that our nonlinear media in this case admit linear magnetic polariza- 
tion. This allows us to reduce the Maxwell's equations to a set of 
amplitude Nonlinear Dirac Equations (NDE). We find two represen- 
tations on NDE- spherical and spinor. In the spherical representation 
we obtain optical votices with spin £ = 1 and in the spinor represen- 
tation, vortices with spin j = 1/2. 
PACS 42.81.Dp;05.45.Yv;42.65.Tg 
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1 Introduction 



The scalar theory of optical vortices was based on the well-known 2D+1 
paraxial Nonlinear Schrodinger equation (NSE) jl], g. A generalization of 
the scalar paraxial theory of optical vortices based on the investigation of co- 
called spatio-temporal evolution equations || f|, [5[ has also been performed. 
The polarization and the vector character of the electric field play an im- 
portant role in dynamics and the stabilization of the vortices. In || this 
problem was discussed for the first time and the existence of polarized vor- 
tices was predicted. To investigate this phenomenon more completely, it is 
necessary to investigated the corresponding vector amplitude equations. In 
this paper, we investigated two type of vector amplitude equations, VNSE 
for three linear othogonal polarized component, and a Dirac representation 
of Nonlinear Maxwell's equations for a combination of linear and circular 
polarized components. 



2 Vector NSE 

The vector version of 3D+1 amplitude equation of electrical field, describing 
the propagation of light in dispersion media with Kerr nonlinearity in a 
dimensionlles coordinate system, moving with group velocity my be writen: 

B A B 2 A ^2 ^ 

2ia— + A ± A-/3— + 7 |A| A = 0, (1) 

where A is the normed amplitude of the electrical field, a = kr^/tov; (3 = 
v 2 k" k; 7 = k 2 rln 2 \A\Q-, are corresponding normed constant, v is group veloc- 
ity, k is carrying wave number, r and t are the spatial and time dimensions 
of amplitude field, k" is the dispersion parameter, n 2 is the nonlinear refrac- 
tife index. We investigate the case (3 = — 1;7 = 1. The constant a has a 
tipical value of a ~ 10 2 — 10 3 in the optical region (a ~ r$k). 

To solve the vector equation (|I|), we use the method of separation of 
variables. The amplitudes of three linear polarized orthogonal vector fields 
A>(x, y, z, t) are represented as: 

A(x,y,z,t)= jA 3 (x,y,z)exp(-iat). (2) 

J=x,y,z 
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In a Cartesian coordinate system, for solutions of the form of (0), the vector 
equation ([]]) is reduced to a scalar system of three nonlinear wave equations: 



a 2 " 1 



2 A l + AA l + l Aj M' = 0; l = x,y,z. (3) 

3=x,y,z 

The system (||]) can be written in spherical variables: 

a 2 A 1 + A r A l + ^Ag !ip A l + ]T (\A^A 1 = 0, (4) 

where I = x,y, z. Next we represent the components of the field as a product 
of a radial and an angular part: 

A l (r } 9, ip) = R(r)Yi(9, <fi);i = x,y, z, (5) 
with the additional constraint on the angular parts: 

\Y x (9,<p)\ 2 + \Y y (6,<p)\ 2 + \Y z (9,<p)\ 2 = const.. (6) 

Multiplying each of the equations by the -^r, and bearing in mind the 
constraint expressed in ([]), we obtain: 

^^(o. + IJ^-^P = /(/+!), (7) 

where i is one number. Thus the following equations for the radial and the 
angular parts of the wave functions are obtained: 

A r R + a 2 R + \R\ 2 R - = (8) 

A eip Y j +£(£ + l) = 0;j = x,y,z. (9) 

As equations shows, the nonlinear term occurs only in the radial 

components of the fields, while for the angular part we have the usual linear 
eigenvalue problem. Going back to set (f|), we see that the separation of 
variables for the spherical functions, which satisfied condition (H), is posible 
only for £ — 1. 



Y x = Ff 1 = sin 9 cos <p\ Y y = = sin 9 sin <p; Y x = Y? = cos 9. (10) 
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By choosing one of these angular components for each of the field components 
we see that the eigenfunctions (|10D are solutions to the angular part of the 
set of equations (13). The nonlinear radial part of equations admit exact "de 



Brogile soliton" solutions [10 



\/2 e iar 

R=^——. (11) 

i r 

The real solutions of the vector amplitude equation (Q) in a fixed basis is 
then: 

A x = $l(R(r)Y x ) = v / 2 Sm( ' ar ^ sin#cos</?cos(at) 

r 

A y = U (R(r)Y y ) = ^/2 sm ^ sin 6 sin (p cos(at) (12) 
A z = ^(R{r)Y z ) = V2 Sm ( aT K osecos(at). 



3 Nonlinear Maxwell's Equations 

We consider the Maxwell's equations for the case a source-free nonlinear 
medium with linear and nonlinear electric polarization and linear magnetic 
polarization (magnetization) f7|, [§]: 

ldB - ldD 

VxB = "9F ;VxB =caF < 13) 

V ■ D = 0;V • B = V • H = (14) 
D = P l + AnP nl (15) 

B = H + AixM h (16) 

where E and H are the electric and magnetic intensity fields, D and B are 
the electric and magnetic induction fields, Pi, P n \ are the linear and nonlin- 
ear polarization of the medium respectively and Mi is the linear magnetic 
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polarization. We apply the slowly varying amplitude approximation to the 
Maxwell equations (|T3|)- (p~6|) : 

E(x,y,z,t)=A(x,y,z,t)e iUot (17) 

H(x,y,z,t) = C(x,y,z,t)e-^ ot , (18) 

where A, C and u>o are the amplitudes of the electric and magnetic fields and 
the optical frequency respectively. After using the Fourier representation 
and condition of slowly-varying amplitudes we obtained the next system of 
Nonlinear Maxwell amplitude Equations(NME), writed in rescaled variables: 



VxA = W2 C-5 ¥ -ft- (19) 



OA d 2 A - - 
V xC = ia 1 A+ — -ifa— + i 1 \A\ 2 A (20) 
at at 2 

V- A = 0;V-(7 = 0, (21) 

where the constants are k\ = we ^ ; k 2 = WA> e ' \i = f^"' ai = ^° r o! Pi = 
k'-r Q /2tl; 7 = r kin 2 \A \ 2 ; 5 = While the nonlinear condition 7 = 1 can 
be satisfied, the constants af, i = 1, 2 have typical values «j ~ 10 2 («j ~ r ki). 
The constants fa have typical values fa ~ 1CT 5 — 1CT 6 << l;fa < Pi and for 
picosecond and sub-picosecond pulses in the transparency region of nonlinear 
optical media they may be neglected. 



4 Dirac representation of NME 

We will neglecting the second derivative in time terms of the eqn.(|TS|)-(^U|). 
as fa « 1. We solve the NME ([19])-(|I]), by the method of separation of 
variables. The slowly varying amplitude vectors of the electric field A and 
the magnetic field C are represented as: 

A(x,y,z,t)=F(x,y,z,t)e iAod (22) 
C(x,y,z,t) = G(x,y,z,t)e iAat . (23) 



5 



Substituting these forms into flTD|)-([2"T]), we obtain the stationary version of 
the NME: 

V x F = -i$ 2 G;V x G = i^F + i-f\F\ 2 F (24) 

V-F = 0;V-G = 0, (25) 

where $1 = ot\ + A«; $2 = SAa — «2 > . When the electric and magnetic 
fields are represented as a sum of a linear and circular polarized component it 
is possible to reduce eqns.(^)- (|25|) to a system of four nonlinear equations. 
Substituting 0: 

! = %F Z - ^ 2 = iF x - F y , ^ 3 = ~G Z ; ^ 4 = -G x - iG y (26) 

into the nonlinear system fl2"4D-(|25D we obtain a stationary nonlinear Dirac 
system (NDE): 

(i-%h + ^ = -^ + ^ 2 ^ < 27 > 
(s + '£)*•- 1*' ^^gi*^*' < 28 > 

— -i— Ur 2 + — ^ = -^3 (29) 
Ox Oy J Oz 

( d d\ d 

— Ui - = -M? 2 * 4 . (30) 

\ox oy J oz 

As may be noted, the optical NDE are significantly different from the NDE 
in field theory. The nonlinear part appears only in the first two coupled 
equations of the system. 

5 Vortex solutions with spin / = 1 

The NDE (f2~7|)-(|3~0l) have both spherical and spinor symmetry only in the 
case where the nonlinear part does not manifest an angular dependence on 

2 

the radial variable \^j\ 2 = F(r). This type of solution may be found 
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using with the following technique: using Pauli matrices, we write the NDE 
system as: 

(ff-P)<P=-t($ 1 +lj2\v j \ 2 )v (31) 

ip ■ P) V = -itirf, (32) 

where a are the Pauli matrices, P — (j^; J^! -§^j are the differential operator 
and rj and (j> are the corresponding spinors. After substituting eqn. ( jjffi) into 
eqn. fl3l]) we obtain: 

2 

(a • • P) V = -M#i + 7 E (33) 

i=i 

When there is no external electric or magnetic field, the operator on the left- 
hand side of eqn. ([33|) is simply the Laplacian operator A. From ( p3[) and 
we obtain: 

2 

$i$2V + $27Y,h\ 2 V + Ar/ = (34) 

3=1 

The scalar variant of these equations has been investigated in many papers 
but exact localized solutions have not been found. In the case / = 1 we look 
for spinors of the form: 

rji = rj(r) cos 6; 7/2 = rj(r) sin 9 exp icf). (35) 

After substituting eqn. (|3~5| ) into equations (0) the following equation de- 
scribing the radial dependence is obtained: 

W+^E l%l 2 ^+ |l + ;f - ptf = °- (36) 

The angular parts are the standard spherical harmonics with / = 1. This 
system has exact vortex c?e Broglie soliton |K| solutions of the form: 



V^exp^V^i^r) , Q7 ^ 
r/(r) = (37) 



if #27 
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6 Hamiltonian representation of the NDE: 
Orbital momentum, first integrals and vor- 
tex solutions with spin j = 1/2 

It is not difficult to show that for the NDE system of eqns. fl2"T|)-(|3"0D a 
Hamiltonian of the form: 

I*/ (38) 

3=1 

may be written. Using this, NDE may be rewritten in the form: 

= etf , (39) 

where e = (— iv\, —ivx, —if 2, — ^2) is the energy operator. Here we investi- 
gate the case where the nonlinear part, as a sum of square module of spinors, 

2 2 

depends only on the radial component = F(r). We also introduce 

here the well known orbital momentum operator L, spin momentum S, and 
the full momentum J, as well as: 

L — fxP;J — L+ -S. (40) 

It straightforward to show that the Hamiltonian of eqn. (|38|) commutes with 
the operators J 2 and J z (the z- projections must be x or y). Using these 
symmetries and the condition that the nonlinearity be centrosymmetric, we 
can solve the NDE equations ( |3"9"| ) by a separation of variables technique. To 
find vortex solution of the NDE ( |39|) with spin j = 1/2 we use the spinor 
representation of these equations. We look for solutions of the form: 

#1 = a(r)ft jim ; ^ 2 = a{r)Sl jlm (41) 

* 3 = ib{r)tl jfm i * 4 = ib(r)Q jl/m (42) 

where fX,/™ is the spherical spinor, 1 + 1=1, and a(r) and b(r) are arbitrary 
radial functions. Using the symmetries of ( |39D and the fact, that the non- 
linear parts depend only on r, we separate variables and obtain the following 
system of equations for the radial part: 
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da(r) 1 + cr.. „ , . . . , . 

Or r 

+ - — -b(r) = iWr) + 7 |a(r)| 2 a(r), (44) 
ar r 

where a = Z(Z + 1) — j(j ; + 1) — 1/4. We find localized solutions of these 
equations only for the angular component corresponding to / = 1 and j = 
1/2. In this case the system (0)"(11) becomes: 

^P- + l a {r) = -$ 2 b{r) (45) 
or r 

^l = § l a( r )+ 1 \a{r)\ 2 a{r). (46) 
As was shown above, this system has exact radial solutions of the form: 

v^exp iy/W^r 
a(r) = — (47) 



b (r) = ~— I 3 + 1 • (48) 



The vortex solutions of the spinor representation have four components and 
this type of solutions is found only for the case where the spin j = 1/2. 



7 Conclusion 

In this paper we investigated two cases of vector amplitude equations of 
electromagnetic field in nonlinear Kerr media. In the first case, for three 
orthogonal linear polarized component of electrical field we obtained exact 
vortex solutions of VNSE with spin 1 = 1. The numerical analysis shows, that 
these solutions are stable at distances comparable to those where localized 
waves with the same amplitudes of the scalar 3D+1 NSE self-focuced rapidly. 
The conditions (3 = — 1 is satisfied only near to Langmuir frequency in plasma 
or near to some of the electron resonances. In the second case we have derived 
a set of Nonlinear Maxwell amplitude Equations (NME) for nonlinear optical 
media with dispersion of the magnetic susceptibility. We show that for the 
case of linear and circularly polarised components of the electric and magnetic 
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fields, the NME is reduced to the Nonlinear Dirac system of equations (NDE). 
Using the method of separation of variables, we have obtained exact vortex 
solutions for this case. We have investigated the NDE in two representations: 
spherical and spinor. In the spherical presentation we obtain optical vortices 
with spin / = 1 and in the spinor representation vortices with spin j = 1/2. 
These vortices are Lorenz invariant, and periodically pulsating with the group 
velocity of solutions of the NME. 
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